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Abstract 

We construct an iterated stochastic integral with fractional Brownian motion with H > 1/2. 
The first integrand is a deterministic function, and each successive integral is with respect to 
an independent fBm. We show that this symmetric stochastic integral is equal to the Malliavin 
divergence integral. By a version of the Fourth Moment theorem of Nualart and Peccati [7], we 
show that a family of such integrals converges in distribution to a scaled Brownian motion. An 
application is an approximation to the windings for a planar fBm, previously studied by Baudoin 
and Nualart [2]. 



1 Introduction 

Let B = {(Bl, . . . , Bf),t > 0} be a multidimensional fractional Brownian motion (fBm) with Hurst 
parameter H > 1/2. We study the asymptotic behavior as k —> oo of multiple stochastic integrals of 
the particular form: 

Y kt := f s-' !H l {1 < Sl< ... <Sq < kt} dB = f f q ■■■ f * s-^dB^ . . . dB«;_\ dB* 

J [0,00)1 Jl Jl Jl 

where t > and each iterated integral is a pathwise symmetric integral in the sense of Russo and 
Vallois [S], and also a divergence integral. Our main result is a central limit theorem for the process 
{Yk*,t > 0}, namely that k converges in distribution as k — > oo to a scaled Brownian motion. Our 
approach uses the techniques of Malliavin calculus, where we express Y k t in terms of the divergence 
integral S, which coincides with the multiple Wiener-Ito stochastic integral in this case. In our proof, 
convergence of finite-dimensional distributions follows from a multi-dimensional version of the Fourth 
Moment Theorem [S] , which gives conditions for weak convergence to a Gaussian random variable 
(see section 2.4). Functional convergence to a Brownian motion is proved by investigating tightness. 
In addition to the proof, we are able to comment on the rate of convergence (which is fairly slow: 
~ (logfc)~2), using a result from Nourdin and Peccati [4] in their recent book on the Stein method. 

The original motivation for this paper was [5] , where Baudoin and Nualart studied a complex- valued 
fBm with H > 1/2. For B = B\ + iBf, B = 1, they studied the integral 

f*dB B f* B\dB\ + B 2 dB 2 . f* B 2 a dB\ - B]dB 2 s 

IT = IR 12 + 1 / R 12 - W 

Jo u s Jo l-t^r Jo \ B s\ 
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When B is written in the form pte l6t , the angle B t is given by the imaginary part of ([T]). For standard 
Brownian motion, a well-known theorem by Spitzer |10j holds that as t — > oo, the random variable 
2#t/(logt) converges in distribution to a Cauchy random variable with parameter 1. We are not aware 
of a corresponding fBm version of Spitzer's theorem. In [2], the functional 

f* BjdBl ~ BldB 2 s 
Z t / — h?TT^ ~ ( 2 ) 



was proposed as an asymptotic approximation for t . It was shown (see Proposition 22 of [5]) that ^f* t 
converges in distribution to a Gaussian random variable, with an expression for variance similar to our 
own result. Their proof also used Malliavin calculus, but did not use the Fourth Moment Theorem. For 
q = 2, since B t = J Q dB s , Z t is asymptotically equal in law to 

, , dBldB\ f* f s dBldB 

and we have a new (and shorter) proof of the result in [2J . 

The organization of this paper is as follows. In Section 2, we give necessary details about the 
theoretical background; this includes a brief discussion of Malliavin calculus, and some remarks about 
integration with respect to fBm, which has been studied extensively elsewhere. The Fourth Moment 
Theorem is also given. In Section 3, we state and prove the main result, which is Theorem 3.2. The 
proof follows in Sections 3.1 through 3.3; then Section 3.4 discusses the rate of convergence. Section 4 
contains a technical lemma which was put at the end due to length. 



2 Notation and Theory 

We use the symbol 1a for the indicator function of a set A. Given a real n-tuple x := (x%, . . . , x n ), we 
can re-arrange the variables in increasing order. We denote this re-ordered n-tuple (scm, . . . , xr n -\). For 
a stochastic process X = {X(a),a £ 1}, we will use the notation X a and X(a) interchangeably. The 
symbol C denotes a generic positive constant, which may change from line to line. 

2.1 Elements of Malliavin calculus 

Following is a brief description of some identities that will be used in the paper. The reader may refer 
to [5] for detailed coverage of this topic. Let X — {X(h), h € H} be an isonormal Gaussian process on 
a probability space (fl,J-,P), and indexed by a Hilbert space H. That is, X is a family of Gaussian 
random variables such that E[-X"(h)] = and E [X(h)X(g)] — (h,g) n for all h,g £ H. 

For integers q > 1, let W® q denote the q th tensor product ofH. We use H &q to denote the symmetric 
tensor product. Given a function / g H® q , we define the symmetrization / g as 

f(x 1 ,...,x q ) = — ^/(av (1)) ...,a; CT(9 )), (3) 

a 

where a includes all permutations of {1, . . . , q}. 

Let {e„, n > 1} be a complete orthormal system in H. For functions f,g£ H Qq and p € {0, . . . , q}, 
we define the p th -order contraction of / and g as that element of 'H® 2 (i~p) given by 

oo 

ti,...,i p =l 

where / ®q g = f (g) g and / ® q g = (/, g) n ® q - While / and g are both symmetric, the contraction may 
not be. We denote its symmetrization by / ® p g. 
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Let H q be the q th Wiener chaos of X, that is, the closed linear subspace of L 2 (S1) generated by the 



r,tk 

vv iciici uiauo ui s\ , unau id, liic uiudcli iiiicai duu 

random variables {H q (X(h)),h £ ||/i||^ = 1}, where H q (x) is the q th Hermite polynomial, defined 

as 



*£ d q 

'dx q ' 

For q > 1, it is known that the map 

I q (h*«)=H g (X(h)) (5) 

provides a linear isometry between the symmetric product space "H 0<? (equipped with the modified norm 
\/ql\\ • an( i where I q (-) is the Wiener-Ito stochastic integral. By convention, Ho = R and 

/o(^) = x. It follows from ([5|) and the properties of the Hermite polynomials that for f,g £ "H 09 we 
have 

E [J, (/)/, (5)] =?!</, 9) H » q . (6) 

Let 5 be the set of all smooth and cylindrical random variables of the form F — g(X(4>i), ■ ■ . , X(cf> n )), 
where n > 1; g : R™ — > R is an infinitely differentiable function with compact support, and <pi G W. 
The Malliavin derivative of F with respect to X is the element of L 2 (fl,H) defined as 

By iteration, for any integer q > 1 we can define the q th derivative D q F, which is an element of 
We let B 9 ' 2 denote the closure of S with respect to the norm || • || D9 ,2 defined as 

\\F\\l,, a =E[F 3 ]+Y i E[\\D i F\\ 2 H9i ]. 

i=l 

We denote by S the Skorohod integral, which is defined as the adjoint of the operator D. A random 
element u € L 2 belongs to the domain of 5, Dom 8, if and only if, 

\E[(DF,u) H }\ < Cu||.F||£a ( n) 

for any F G D 1,2 , where c u is a constant which depends only on u. If u G Dom 5, then the random 
variable 6(u) G L 2 (J7) is defined for all F 6 D 1 ' 2 by the duality relationship, 

E[FS(u)} = E[{DF,u) n ] . 

This is sometimes called the Malliavin integration by parts formula. We iteratively define the multiple 
Skorohod integral for q > 1 as 5(5 q ^ 1 (u)), with 8°(u) — u. For this definition we have, 



E[F5 q (u)} =E[(D q F,u) 



where u E Dom 5 q and F € D 9 ' 2 . The adjoint operator 5 9 is an integral in the sense that for a 
(non-random) h £ H Qq , we have 8 q (h) = I q (h). 

We will use the following hypercontractivity property of iterated integrals (see [7], Theorem 2.7.2, 
or [5], Sec. 1.4.3 for complete details). Let / £ T-L 0q and p > 2. Then there exists a positive constant 
C Pi q < oo, depending only on p and g, such that 

E[|/,(/)n<C M (E[/ 9 (/) 2 ]) f . (7) 
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2.2 Fractional Brownian motion 

Fix T > and an integer d > 1. Let B = {B u < t < T} = (B\, . . . , Bf) be a rf-dimensional fBm, that 
is, each B\ is an independent, centered Gaussian process with Bq — and covariance 

E [BlBl\ := R(s, t) = \ (s 2H + t 2H -\s- t\ 2H ) 

for t, s > 0. We assume that | < if < 1. We will use the following elementary properties of R(s, t): 
(R.l) R(s, t) = R(t, s); and for any e > 0, R(s + e, t) > R(s, t). 
(R.2) There are constants 1 < Co < c\ < 2 such that Co(st) < R(s,t) < c\(st) H . 
(R.3) As an alternate bound, if s < t then the Mean Value Theorem implies 

R(s, t) < s 2H + t 2H - (t - s) 2H < s 2H + st™- 1 . 

Let £ denote the set of K— valued step functions on [0, T] X {1, . . . , d}. Note that any / = f(t, i) € £ 
may be written as a linear combination of elementary functions e\ — l[o,t]x{fc}- Let Sjd be the Hilbert 
space defined as the closure of £ with respect to the inner product 

k j\ _ wtnk Rii 



elel) =E[B«Bi]=R(s,t)6 kj , 
i fid 

where Skj is the Kronecker delta. The mapping e\ i— > B k (t) can be extended to a linear isometry between 
S)d and the Gaussian space spanned by B. In this way, {B(h), h £ fid} is an isonormal Gaussian process 
as in Section 2.1. 

Let aH = H(2H — 1). It is well known that we can write 

R(s,t) = a H [ f \T]-e\ 2H - 2 d v d6. (8) 
Jo Jo 

Consequently, for f,g G £ we can write 

(/, g)* d = ol H y; [ T [ T f( s , i) g (t, i)\t ~ s \ 2H - 2 d S dt. (9) 

i=1 Ja Jo 

We recall (see [5], Sec. 5.1.3) that F)d contains the linear subspace of measurable, ]R-valued functions 
ip on [0, T] X {1, . . . , d} such that 



Y] [ ( \<p(e, i)\ \<p{t, i)\\t- s\ 2H - 2 ds dt < oo. 
. 1 J o Jo 



o Jo 

We denote this space \Sjd\- Let |i^' s | be the space of symmetric functions / : ([0, T] x {1, . . . , d}) q 
such that 

d . 

V / \f(( Vu h),...,( Vq J q ))\\f((0 u h),...,(O q ,i q ))\\v-O\ 2H - 2 dr)de <oo. 
Then \S) d ' s \ c fi &q : an d for /,j6 |£>J S | we can write (@| as 

p 



f® P 9 = a p H J2 [ /((r 7 ,fc),(t 1 ,i 1 )) ff ((0,fc),(t 2 ,i 2 ))f[|r /j -^| 2H - 2 ^^, 



(10) 

where 

(rj,k) = (r)i,k),...,(r) p ,k); (8, k) — [0i,k), . . . , [6 P , k)\ (ti, ii) = (ii, h), . . . , (t q - p ,i q - p ); and 

(t2, b) = (*g-p+l, *g-p+l)i • ■ ■ : (*2(g-p)j*2(g-p))- 
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2.3 Stochastic integration with respect to fBm 

Let F = g(B(<fii), . . . , B(cj) n )), where n > 1, g : M. n — > R, fa £ fid, and g is a smooth function as in 
Section 2.1. The Malliavin derivative of F is an element of S)d (which is isomorphic to the product 
space (Sji) d ), and we can write D — (D^\ . . . , D^ d >), where 

n ft 

where we use the notation F — D^F{t). We define the 'component integral' #W as the adjoint of 
DW, and use the notation 

<S w (u) = / UtSBl; and (11) 



«(«)= f*=r^(i 

^0 i=l 



where u € Dom <jW c L 2 (0,i5i) for every i = 1,. . . ,d implies u £ Dom S C L 2 (f2, fjd). 

The pathwise stochastic integral with respect to fBm with H > 1/2 has been studied extensively 
[H G3 [5]. For our purposes, we will use the symmetric Stratonovich integral discussed by Russo and 
Vallois [9]: 

Definition 2.1. For some T > 0, let u — {u t ,0 < t < T} be a stochastic process with integrable 
trajectories. The symmetric integral with respect to the fBm B is defined as 



/ u s dB s = lim — / u s (B {s+e)m - £ (s _ e ) V0 ) ds, 
Jo do 2e J 



where the limit exists in probability. 

This theorem was first proved in pQ. 
Theorem 2.2. Let u — {u t ,t > 0} be a stochastic process in B 1,2 (f)i) such that, for some T > 0, 

cT i-T 



E 



E 



< oo; 



/ / \u t \ \u s \ \t - s\ 2H ~ 2 ds dt 
Jo Jo 

[ \D t u g \ \D s u v \ \t - s\ 2H - 2 \6 - ii\ 2H - 2 du dt d6 drj 



< oo; 



and 



T r T 



o Jo 



\D s u t \ \t - s\ 2H 2 ds dt < oo a.s. 



Then the limit of definition 2.1 exists in probability, and we have 



u t dB t = \ u t SB t + a H 



T r T 



D s ut\t- s\ iU ~ i da dt 



where an — H(2H — 1). 



G 



2.4 The Fourth Moment Theorem 

Theorem 2.3. Fix integers n > 2 and d > 1. Lei | (fi^i ■ ■ ■ , /j*^ > l| ^ e a sequence of vectors 
such that G H ™ /or eac/i fc and i = 1, . . . , d; and 



lim niy/j ; ||l/®n = lim 

A:— foo /c— foo 



4 /, 



P (fc) 



L 2 (£l) 



= Q», Vi = 1, . . .,d; 



lim E 

k—too 



dj, VI < i < j < d. 



Then the following are equivalent: 



(i) As k —¥ oo, the vector (l n (f^), . . . , I n {f^)^j converges in distribution to a d- dimensional Gaus- 
sian vector with distribution A/"(0, C4); 

(ii) For each i = l,...,d, Inifi^) converges in distribution to Ni, where Nj, is a centered Gaussian 
random variable with variance Cu; 



(Hi) For each i = 1, ... ,d, 



lim 1 

k— f 00 



M/, 



f (fe) 



= 3Ci 



(iv) For each i = 1, . . . , d, and each integer 1 < p < n — 1, Iinifc- 



(fc) 



(fc) 



^®2(n-p) 



0. 



This first version (which was 1-dimensional) of this theorem was proved in [7]. Since then, other 
equivalent conditions have been added [4,6. The multi-dimensional version stated above was proved 
by Peccati and Tudor [8]. A key advantage of this theorem is that, unlike the standard method of 
moments, it is not necessary to know about moments of any order higher than four. 



3 Main result 

Fix q > 2. For t > and integer k > 2, define 

Y kt ^ qH dB l sl . . . dB«;_\ dB« , 

where the stochastic integrals are iterated symmetric integrals in the sense of Definition 2.1. Theorem 
2.2 and the diagonal structure of Y^t allow us to identify the pathwise and Skorohod integrals. 

Lemma 3.1. For each q > 2, we have 

Y„ = j* £■■■ jT 8-« B 8B^ . . . 5Bj;_\ SBl q . (12) 

Proof. This follows from iterated application of Theorem 2.2, where the correction term is zero due to 
independence. Indeed, in the notation of (ITT1) . this is 

Y kt (s-« H l { i< Sl< ... <Sg < kt} ) . 

□ 

Following is the main result of this section. 
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Theorem 3.2. For t > 0, define 

X k (0) = 0; X k (t) = -jk=, t > 0. 

Then as k — > oo, the family {X k (t),t > 0} converges in distribution to the process X = {X(t),t > 0}, 
where X is a scaled Brownian motion with variance o~ 2 t, and 

o-\ = a H \ x~ 2H R(l,x){l- x) 2H - 2 dx; and for q > 2, (13) 
Jo 

a\ = a^ 1 / x~ qH {l - x q ) 2H - 2 / R(x 2 ,y 2 ) TT \ Xl - y t \ 2H ~ 2 dx 2 dy 2 . ..dy q - X dx g , (14) 
Jo Jm i=2 

where M — {0 < x 2 < ■ ■ ■ < x q ; < y 2 < ■ ■ ■ < y q -\ < 1}. 

The proof of Theorem 3.2 follows the lemmas in Sections 3.1 and 3.2. Our first task is to investigate 
the covariance (Section 3.1), then verify two other conditions for weak convergence (Section 3.2). 



3.1 Convergence of the covariance function 

Let A = {1 < si < ■ ■ ■ < s q < k f }, and B = {(ii, . . . , i q ) = (1, . . . , q)}. Lemma 3.1 allows us to write 
Y k t = 5 q (f k t), where f k t : ([0, oo) x {1, . . . , q}) q -> R is given by 

fk* ((Sl,il), (Sq,iq)) = Sq qH l A (si, ■ . .,S q )l B (il, ■ ■ .,«,). (15) 

Here, f k t £ S^ q , where fj := S) q is the Hilbert space associated with a q— dimensional fBm (see Section 
2.2). Clearly, f k t is not symmetric. Instead, we will work with the symmetrization defined in (J3]): 

fk' = ^^s~{>" l A (s a{1) , . . . , s a{q) )l B (i a (i), ■ ■ ■ ,i*( q )), (16) 

where a covers all permutations of {1, ... , q}. This gives equivalent results, by the relation I q (f) = I q (f) 
(see [5], Sec. 1.1.2). 

By definition f k t is nonzero only if 1 < So-fi) < • • • < s a(q) < k* and (V(i)i ■ ■■■> = (1> • • • > q)i 

hence it is possible to express f k t without a sum. Let a be an arbitrary permutation of {1, . . . , q}, and 
let A a = {1 < s CT (i), < • • • < s CT ( 9 ) < k }. Since the sets {^4o-} form an almost-everywhere partition of 
[1, A:*] 9 , we can write (fT6|) as 

A* = ~\-\q) Hl M {(Sl,il),---,(s q ,iq)), (17) 

where st q \ — max{si, . . . , s q }, and the set A\ is defined by the following condition: when s\, . . . , s g are 
arranged in [1, fc*] such that S(i) < • • • < S( 9 ), then (i(i), • • • , i( q )) = (1, . . . , q). 

In the next three results, we check the conditions of Theorem 2.3 for 5 q (f k t). 

Lemma 3.3. For each q > 2 and t > 0, 

to* = hm E[X k {t) 2 ] 
exists, where a 2 is given by (|13[) and (I14p for q = 2 and q > 2, respectively. 
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Proof. Since is deterministic, we use ((6]) and (|9|): 



g! lo : 



/ (r (?)S(g) )-^l Al (r,i)l Al (s,j)TT|r,- S ,| 2 ^ 2 dr ds, (18) 



where (r, i) = ((ri, . . . , (r 9 , and similar for (s, j). To evaluate (ITBl . we decompose [1, k t ] 2q into 
the union of the sets {A a x A a >}, which form a partition almost everywhere. Since 1a 1 (r, i) is nonzero 
only if r CT m < • • • < and (vm, • ■ • j V(g)) = (1> - • • > 9)1 an d similar for l^ 1 (s,j), it follows that we 
integrate only over the diagonal sets, that is, when a — a' . Hence, (fT5)l can be integrated as a sum of 
q\ equal terms, and we have 



Ti — Si I dri ds\ . . . dr q ds q , (19) 



where the integral is over the set 

A = {1 < n < ■ ■ ■ < r q < k\ 1 < si < • • • < s q < fe*} . 
Integrating over r\, s\, we have by L'Hopital, 

9-1 r r 9-1 



q —± /> p H 

lim / {r q s q )- qH \r q - s q \ 2H - 2 / i?(r 2 ,s 2 ) TT |r< - S;| 2iy - 2 dr 2 ds 2 . ..dr q ds q 

:^oo logfc 7[l,fe*]2 y^t' f = 2 

= lim ifc'a^ 1 / (r^r^fe* - r,) 2 ^ 2 / R(r 2 ,s 2 ) TT |r ( - s 4 | 2 ^ 2 <*r 2 da a • • • dr 9 , 



i=2 

where the set .A' = {1 < r% < ■ ■ ■ < r q , 1 < s 2 < ■ ■ ■ < s g -i < k 1 } (A' is empty if q = 2). Using the 
change of variable = a;,-, Sj = y,, this may be written 

r 1 r 9-1 

lim ta"- 1 / 1-^(1 - x q ) 2H - 2 / i?(x 2 ,y 2 ) TT \ Xi - Vl \ 2H ~ 2 dx 2 dy 2 . ..dy q _ x dx q (20) 

fe^oo Ji J M fj- 



i=2 

-1 f 1 f TT 

ia^ 1 / x~ qH (l - x q ) 2H ~ 2 j R(x 2 ,y 2 ) TT \x t - yi\ 2H ~ 2 dx 2 dy 2 . . . dy q -i dx q , 
Jo Jm i=2 



where M is as in (fill) for q > 2, and we have (fT^ if g = 2. To show (fT3|) and (fTl|) are convergent, we 
use properties (R.l) and (R.2), so that 



a 2 = a H f x- 2H {l - x) 2H - 2 R(l, x) dx < Cl a H f x~ H (l - x) 2H - 2 dx 
Jo Jo 



< 00 



and for q > 2 

o- 2 < a H [ x- qH {l- x^^Ril^XqY^dxq < c\~ r a B [ x- H (1 - x q ) 2H - 2 dx q < 00. 
Jo Jo 

This concludes the proof. □ 
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Lemma 3.4. Let < r < t. For each q > 2, 

lim E[X k (t)X k (T)] =o\t- 

and consequently lim; i: _>. 00 E [X(s)X(t)} = <J 2 (s A t) for all < s, t < oo. 
Proof. 

E [X fe (t)X fe (r)] = E [(X k (t) - X fc (T) + X fe (r))X fe (r)] 
1 



logfc 



E[(F fet -r fc .)^]+E[X fe (r) 2 ] , 



where E [JTfc(r) 2 ] — > <t 2 t by Lemma 3.3. Note that Y^t — Y^ = <5 9 (/fct) — 5 q (fk T ), where, recalling the 
notation of (fTTj) . 

<5 ? (A*) - S«(frr) = [ * l Al (s,i)5B s - [ *i Al fai) S B a 

' 5b£1 ■ ■ ■ SBlr 1 ! SBl*\ - / ■ ■ • / -4tf*I#> ■ ■ ■ SBir 1 ! *B$ 



I Off S(l) s (9) / / / S(l) 8(3-1) S(<r) 

1 9.S (?) Jl Jl Jl q.S {q) 



fc* /■»(„) rs {2 ) i 



• • -(5B (9_:L) 



Hence, we can write Y^t — Y^t — (5 9 (/a/c)i where 



/a/c — _ff 1 -4i 1 {fc T <S(,)<fe t } — /* t l{fc T <«( 9 )<fc t }- (21) 



With this notation, it follows that 



1 -E [(Y fc * - Y fe r)Y feT ] = (/ Afe , /*, 



log k log fc \ ' / 

i r q 
= / (ww)"'* Wr, iJlA^s, j)l{i<. w <fc-<r (<) <fc*} II ^ " s e \ 2H ~ 2 ds dr. 

As in Lemma 3.3, we decompose [1, k t ] 2q into the union of the sets {A a x Since (r, i) is nonzero 

only if r a m < ■ ■ ■ < r a t g ) and (i CT (i), • • • , Wq)) = (lj • • • j <z)i an d similar for l Al (s, j), it follows that we 
integrate only over the diagonal sets, that is, when o = & ' . Hence, we have q\ equal terms of the form 



k* pk T 



f q 

„,,,,./ i / hs,) _,ff l{ ri <...< r jl{ Sl <...< S! ) [[ |r« - S< r^(ls dr. (22) 

q\ log fc J k r J x J[ l fc t]2,-2 ^ 

By (R.l) and (R.2), for each rg < r g , si < s q , we have the estimate 

an J J V{i-X) - S{i-\)\ ~ dry-i) dsy-i) 

<a H / \r-s\ 2H - 2 dr ds = R(r {q) ,s (q) )<c 1 {r qSq ) H . (23) 
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It follows that 

} ±-E[(Y kt -Y k r)Y k ,] < f [ k (r q s q )-« H R{r q ,s q y- l \r q - s q \ 2H - 2 dr q ds q 
log k log k J kT J 1 

q rk* r k T 



~ loifc J k r l (r q s q )- 2H R(r q , Sq )\r q ~ s q \ 2H - 2 dr q ds q . 



Using the change-of-variable s q — k T x, r q — k T y, this is bounded by 

(xy)- 2H R{x, y)(y - x) 2H - 2 dx dy. 



c ^ - 1 



log k J t 

Using (R.3), we obtain the estimate 



ii 



c r**~ T ,1 



log k J 1 

where 



{y- 2H (y - x) 2H - 2 + x l - 2H y-\y - x) 2R - 2 ) dx dy, 



fe /-i r 2 rv r k 

2H/„. „\2H-2j„,.^ / „,-2H / / „\2H-2 > , , , / /„, -,\-2 



V "':// -rr" 2 ,l.r,l, r J V~™ J (y xf» 2 dx dy • J ( „ - I) 

/2 poo 
y~ l dy + C J y~ 2 dy < oo, 



and 



rl y- x x x - 2H (y - x) 2H - 2 dx dy < f y- 1 f (y - xf H ~ 2 dx dy + /* (y - l) 2H - 3 dy 

./ii Jl Jo J 2 

.2H-3, 



<C^ y- 1 dy + J 



y A dy < oo. 

Hence, this term vanishes and Lemma 3.4 is proved. □ 
3.2 Conditions for weak convergence of {X k (t)} 

In the next two lemmas we verify additional properties of {X k (t)}. In Lemma 3.5 we check condition 
(iv) of Theorem 2.3, and Lemma 3.6 is a tightness result. 

Lemma 3.5. Fix q > 2 and t > 0. For each integer 1 < p < q — 1, 

lim (logfc)~ 2 ||/ fet ®p /fet|| 2 ,«2(,- P ) = 0. 

Proof. Let 1 < p < — 1. To compute the p th contraction of f k t, we use (fT0|) . 

fkt®pfkt = T-%2 / ( r (?) s (g)) _9Hl A 1 (r,i)lA 1 (s,j)TT|^-s £ | 2H - 2 dri d Sl ...dr p ds p . (24) 

J[l,fc*p ^ 

Using (|24l) . we want to compute 

ll/fe* ®p /fc* || fi®2(,-p) = </fe* ®p /fe*)/fe* ®p /fe* ) 
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2q r „w P 

<*H ( , , \~ qH \ 4 T T / 1 ,, , m2i?-2 



4 



(?!) 



[l,fc*] 4 9 



( r (9) s (9) r ( 9 ) s (<;)) (iAi) 4 !! (ki - S *IK - s ilf 



t=l 



II {V l -r' l \\s l -s' l \f H ~ 2 dr dsdr'ds'. (25) 



As in the proof of Lemma 3.3, we view integration over the set [l,fc*] 49 as a sum of integrals over 
various cases corresponding to the orderings of the real variables r±,...,r q (as in Lemma 3.3, the 
variables s,r',s' must follow the same ordering). Up to permutation of indices, each integral term has 
the form 

A /(n 3 )%)^) s ( 8 ))"^n(^- s *ll^-^l) 2H " 2 II (Iri-rillsi-s^-'drdsdr' ds', 

Jq i=l i=p+l 

(26) 

where Q = |l < < ■ ■ ■ < < fc*; . . . ; 1 < s'^ < ■ ■ ■ < s'^ < fc* |. To evaluate (l2l)l) . there are two 
cases to consider. The first case is if r^ e {rj., . . . , r p }, that is, (j2"o]l contains the terms |f( 9 ) — S( g ) |, |r| ? s — 
In this case, using (|2"3")l we can bound (121)1) by 



(«0 



( r (9) s (?) r (9) s '(g)) 9 ( fl ( r fe)> s (9))^( r W> S (?))) (- R ( r (9)' r (?)) i? ( s (9)' s '(g)) 
[l,fc*] 4 v y v 



/ \ 2H-2 

X (l r (9) " S (?)l K,) - S (9)l 



dr (g) ds (g) dr' (g) ds' (g) 



< C / (rs/s')" -R(r, r') i?(s, s') (|r - s| \r' - s'|) dr ds dr' ds', (27) 

where we used (R.2) in the last estimate. The second case is the complement, that is, r/q\ G {r p+ \ 1 . . . , r q \ 1 
so that PS]) contains the terms \r^ — r',-.\, |s( 9 ) — siA. If this is the case, then (|26p is bounded by 



4 



(?0 



[i,fe*] 



( r (?) S (?) r ( g ) S (9)) ( i? ( r (9)' S (9)) i? ( r '( 9 )' S '( 9 ))) P ( i? ( r (9)' r ( g ))- R ( S («)' S ( 9 ))" 

/ \ 2H-2 

X (l r (9) - S (l)\ \ r U - S {q)\) dr (l) dS (l) dr '(q) dS '(q) 



<C {rsr's')~ 2H R(r, s) R(r',s') (|r - r'\ \s - s'\f H ' 2 dr ds dr' ds' . (28) 

J[i,k t ] i 

The result then follows by a change of variable and applying Lemma 4.1 to (|27|) and (|28|) . □ 

Lemma 3.6. There is a constant < C < 00 such that for each k > 2 and any < t < t < 00 we have 

E[\X k (t)-X k (r)\ 4 ] <C{t-r) 2 . 
Proof. Based on the hypercontractivity property (J7J), it is enough to show 

E[\X k (t)-X k (r)\ 2 ] <C(t-r). 
Using the notation of (|2"Tj) . we can write 

E [\X k (t) - X k (r)\ 2 ] = J-^E [\Y kt - n, \i = A (; Afc , ; Afc )^ 



TTT / { r {q) s (q)Y qh Wr,i)lAi(s,j)l { ^< < fct} []|r £ -s^| 2H - 2 dsdr. 
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In the same manner as (|22[) . this can be decomposed into a sum of q\ equal terms of the form 



q\ log fc 

Similar to Lemma 3.5, we use (|23j) and a change-of-variable to obtain 



i /.1 

logfc 



< At I I (xy)- 2H R(x, y)\x - y\ 2H - 2 dx dy 



Without loss of generality, assume x < y. By (R.3), we have the estimate 

( 1 (xy)- 2H R(x,y)\x-y\ 2H - 2 dxd y =-^- [ f (xy)- 2H R(x,y)\x - y\ 2H - 2 dx dy 



lOgfc J k r- 1 ./,, - ■ l0U/.- 

C 



< r^j I 1 f ?T 2ff (» - ^) 2H - 2 + x'- 2H y-\y - x) 2H ~ 2 dx dy 
log A; J k r-t Jo 

<At t y-^y^-' + y-'dyKCit-r). 



logfc J kT - 

This concludes the proof. □ 

3.3 Proof of Theorem 3.2 

Fix integers q > 2 and d > 1, and choose a set of times < t\ < ■ ■ ■ < td- Lemmas 3.3 and 3.4 
show that the random vector sequence {(Xk(ti), . . . ,X)~{td)) , k > 1} meets the covariance conditions 
of Theorem 2.3. Moreover, Lemma 3.5 verifies condition (iv) of Theorem 2.3. Therefore, we conclude 
that as k — > oo, 

(Xfc(ti), • • • , X k {t d )) -A {X{h), X(t d )) , (29) 

where each X(ti) has distribution 7V(0,ct 2 ^), and E [X(U)X(tk)} = <r 2 (U A tk) for all 1 < i, k < d. By 
Lemma 3.6, the sequence {Xk(t)} is tight, hence it follows from (|2"9")l that the sequence converges in 
the sense of finite-dimensional distributions, and we conclude that the family {Xk{t),t > 0} converges 

in distribution to the process {X(t),t > 0} = {a q Wt,t > 0}, where Wt is a standard Brownian motion. 
This concludes the proof of Theorem 3.2. 

3.4 Rate of convergence 

Let t > be fixed. By Theorem 3.2, it follows that the sequence {Xk(t), k > 1} converges in distribution 
to a random variable N(t), where N(t) ~ JV(0,a 2 t). Recent work by Nourdin and Peccati [4] has 
produced a stronger form of the Fourth Moment Theorem for the 1-dimensional case, that is, that the 
conditions of the Fourth Moment Theorem also imply convergence in the sense of total variation (as 
well as other metrics - see Theorem 5.2.6). The result below follows from Corollary 5.2.10 of [3]. 

Proposition 3.7. Let t > 0. Then for sufficiently large k, there is a constant < C < oo such that 

d T v(X k {t),N{t))< ( 



%/logfc 

where drv(', ■) is total variation distance. Hence Xk(t) converges as k — » oo to Gaussian in the sense 
of total variation. 
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Proof. The result follows from an estimate in |2] (Cor. 5.2.10): 



try (XM. iVW) < 2. E ^ - 3E ^ + 2 - ffi . (30) 

To simplify notation, we will assume t = 1. To help interpret this estimate, the following identity is 
computed in [1] (see Lemma 5.2.4): 

E [X k (l) 4 } -3E [X fe (l) 2 ] 2 = -^-^^(p!)2Q 4 ( 2(7 - 2p)!||/ fc ® p /*||* 8a(4 _,>. (31) 

From Lemma 3.5, we know (log k)~ 2 \\fk ®p /fe||^®2(,- P ) — > at a rate C/logfc, hence it follows the first 
term of (|3U|) is of order C(logfc)~2 . The second term depends on the covergence rate of (|T5| . In the 
proof of Lemma 3.3, convergence follows from a limit of the form E [Y^l/logfe. By L'Hopital's rule, it 
follows the rate of convergence has the form C / log k, hence the first term controls. □ 



4 A technical lemma 

Lemma 4.1. FixT > 0. Let 1/2 < H < 1, and for nonnegative x, y, let R(x, y) = | (x 2H + y 2H — \x — y\ 2H ) ■ 
Then there is a constant < K < oo such that 

(xyuv)~ 2H R(x,y) R(u,v)\x - u\ 2H ~ 2 \y - v\ 2H ~ 2 dx dy du dv < KlogT. 

i] 4 

Proof. In the following computations, we will obtain estimates based on the order of integration. Due to 
the symmetries of the integral, it is enough to consider four distinct cases. We will make frequent use of 
(R.3), and for a second estimate, note that for x < y < u we can write (u — x) 2H ~ 2 < (u — y)~ a (y — x)~P, 
where a, (3 > satisfy a + f3 = 2 - 2H. 

Case 1: x < y < u < v We can write 

1 rv ru ry 

(uv)- 2H R(u, v) / y~ 2H {v - y) 2H ~ 2 / x- 2H R(x, y)(u - x) 2H - 2 dx dy du dv 



i_ J x 

T T 



<C f f\uv)- 2H R(u,v) [" y- 2H (v-y) 2H - 2 (u-y)- a f\y - x)~^ + x 1 - 2H y 2H - 1 {y - x)~ p dx . . . dv 
<C I / {uv)~ 2H R(u, v)(v - u)- a I y 1 - 2 "-?^ ~ y)-' 3 - a dy du dv 

f 1 pv 



<C (uv)- 2H R(u, v)(v - u)- Q u- /3 du dv 

< C / tT 2ff f u- 2ff (v - u)- a (u 2H + uv 2 "- 1 ) du dv 

<C J v~ l dv < KlogT. 
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Case 2: x < y < v < u For this case, we use constants a,0 > such that a + f3 = 2H — 2, and 
7, S > such that 7 + 5 = a. 

1 [-U r-u r-y 

{uv)- 2H R{u, v) I y- 2H (v - y) 2H - 2 J x' 2 " R(x, y)(u - x) 2H - 2 dx dy dv du 



1_ } _1_ 

T T 



< C j j\uv)- 2H R{u, v) ^ y 1 - 2 "-^^ - y) 2H - 2 {u - y)-° dy dv du 

< C [ [ (uv)- 2H R(u, v)(u - v)-^ [ y 1 - 2 "-?^ - y) 2H - 2 - 5 dy dv du 

J _L_ J 1_ J 1_ 

< C I' u- 2H j" V - 2H -P- S (v 2H + vu 211 - 1 ) (u - v)-^ dv du 
<C J i u- 1 < K\ogT. 

Case 3:x<u<y<v 

fl f ' (yv)- 2H (v - y) 2H - 2 ( V u- 2H R{u,v) x~ 2H R(x,y)(u - x) 2H ~ 2 dx du dy dv 

i_ I _L / j_ / j_ 

T T T T 

< C C f\yv)- 2H {v - y) 2H - 2 I" u- 2H (u 2H + m; 2 "" 1 ) (u 2 "- 1 + y 2 "" 1 ) du dy dv 

< C f (\yv)- 2H (v - y) 2H - 2 (u 2 *- 1 + y 2 "" 1 + v 211 ' 1 + u l ~ 2H (vy) 2 "- 1 ) du dy dv 

< c f f M" 2 > y) 2H - 2 (y 2H + w 2 "- 1 ) d y dv 

j — j — 

< C [ v- 1 dv < K log T. 



Case 4:x<v<u<y 
1 rv 



l\uy)- 2H I" v- 2H R(u, v)(y - v) 2H - 2 f x' 211 R(x, y)(u - x) 2H - 2 dx dv du dy 

1 / J_ I ± / J_ 

T T T T 

< C j f\uy)- 2H v- 2H R(u,v)(y-v) 2H - 2 (u-v)- a f x~ 2H (x 2H + xy 2 "' 1 ) (v - x)- fj dx dv du dy 
<C I I {uy)- 2H (y - u)~ a I v- 2H (v 2H + vu 2 "- 1 ) (u - v)-*-? (v 1 ^ + v 2 - 2H ^ y 2H - 1 ) dv du dy 



1 r y 



<C I I (uyy^iy-u)-*^ 21 *- 13 +y 2H - 1 u 1 - f! ) du dy 

J — J — 

<C [y-^iy'-^ + y 2 "- 1 ) dt 
<C J y- 1 dy< K log T. 



□ 
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